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Magneto-elastic oscillations and the damping of crustal shear modes 
in magnetars 
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ABSTRACT 

In a realistic model of magneto-elastic oscillations in magnetars, we find that crustal shear 
oscillations, often invoked as an explanation of quasi-periodic oscillations (QPOs) seen af- 
ter giant flares in soft gamma-ray repeaters (SGRs), are damped by resonant absorption on 
timescales of at most 0.2s, for a lower limit on the dipole magnetic field strength of 5 x 10^'^G. 
At higher magnetic field strengths (typical in magnetars) the damping timescale is even 
shorter, as anticipated by earlier toy-models . We have investigated a range of equations of state 
and masses and if magnetars are dominated by a dipole magnetic field, our findings exclude 
torsional shear oscillations of the crust from explaining the observed low-frequency QPOs. In 
contrast, we find that the Alfven QPO model is a viable explanation of observed QPOs, if the 
dipole magnetic field strength exceeds a minimum strength of about several times lO^^G to 
lO^^G. Then, Alfven QPOs are no longer confined to the fluid core, but completely dominate 
in the crust region and have a maximum amplitude at the surface of the star 
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1 INTRODUCTION 

The observation of giant flares in soft gamma-ray repeaters (SGRs; 
compact objects with very strong magnetic fields or magnetars 
( [Duncan & Thompson] 1992) ) may open a gateway towards the ex- 
citing field of neutron star seismology. In the decaying X-ray tail 
of two such events, SGR 1900+14 and SGR 1806-20, a number 
of long-lasting, quasi-periodic oscillations (QPOs) have been ob- 
served (see |Israel et al.H2005^ and |Watts & Strohmayer| ( |2007| l for 
recent reviews). Early models interpreted the observed QPO fre- 
quencies as directly related to torsional shear oscillations of the 
solid crust of a neutron star excited during a giant flare event (see 
Duncan] ((T998l);|Strohmayer & Watts] p005l>;[PiroH2005) ; |Sotani| 
et al.] POO?^ ; Samuelsson & Andersson] ( [2007^ , and references 



therein) raising hopes that through their identification, the physi- 
cal properties of the crust could be probed (]Steiner & Watts 2009 1. 
Due to the extremely strong magnetic fields present in magnetars, 
however, a self-consistent model that includes global magnetohy- 
drodynamic (MHD) oscillations interacting with the shear oscilla- 
tions of the crust is required dLevin' (2006), Glampedakis & An-] 
]dersson] pOOet , ]Levin] p007|>, [Lee ( 2007 2008) ). In a highly sim- 



plified model, ]Levin (20071 showed that shear oscillations will be 
absorbed by an MHD continuum of Alfven oscillations, while long- 
lived QPOs may still appear at the turning points or edges of the 
continuum. 



Sotani et al.] J2008b| and subsequently ]Cerda-Duran et aT] 
( ]2009^ (see also Colaiud a et al.] ( 12009] l), using a more realis- 
tic, general-relativistic MHD model but still ignoring an extended 
crust, found two families of Alfven QPOs related to turning points 
of the frequency of torsional Alfven waves near the magnetic pole 
and inside a region of closed magnetic field lines near the equa- 
tor Each QPO family consists of two sub-families differing by 
their symmetry behaviour with respect to the equatorial plane. The 
results of the numerical simulations were explained by a semi- 
analytic model based on standing waves in the short-wavelength 
limit l lCerda-Duran et al.][2009:) . The Alfven QPO inodel is very 
attractive, because it reproduces the near-integer-ratios of the ob- 
served 30, 92 and 1 50 Hz frequencies in SGR 1 806-20, at magnetic 
field strengths expected for magnetars. In this model, the observed 
SGR QPOs, rather than probing the crust, yield information on the 
magnetic field and the compactness of the star. 

The omission of an extended crust in the previous studies of 
]Sotani et al.]p008bt , |Cerda-Duran et al.| ( |2009l ), and ]Colaiuda"etar] 
P009| l can be considered as a limiting case of a very strong mag- 
netic field. For intermediate magnetic field strengths, however, an 
understanding of magnetar oscillations requires the inclusion of 
crust-core coupling. In this Letter, we present the first such real- 
istic simulations of coupled, magneto-elastic oscillations. We use a 
general-relativistic framework, a dipolar magnetic field, and a tabu- 



2 Michael Gabler, Pablo Cerdd-Durdn, Jose A. Font, Ewald Miiller and Nikolaos Stergioulas 



lated equation of state (EOS) for dense matter. The numerical sim- 
ulations are based on state-of-the-art Riemann solver methods for 
both the interior MHD fluid and the crust. A recent study by van 
|Hoven & Levin| l [2010 ) also takes entanglement of magnetic field 
lines into account, thereby generalising the toy model of |Levin| 
( |2007) . Some first results on coupled crust-core oscillations also 
appeared in |Kokkotas et al.|(2010[ >. 

We use units where c — G = 1 with c and G being the speed 
of light and the gravitational constant, respectively. Latin (Greek) 
indices run from 1 to 3 (0 to 3). 

2 THEORETICAL FRAMEWORK 

The present study of torsional oscillations of magnetars is based on 
a numerical integration of the general relativistic MHD equations. 
As in [Cerda-Duran et al. (2009 ), who considered purely Alfven 
oscillations of the fluid core, assume (i) a zero temperature EOS, 
(ii) axisymmetry, (iii) a purely poloidal magnetic field configura- 
tion, (iv) the Cowling approximation, (v) a spherically symmetric 
background, and (vi) small amplitude oscillations. Because of as- 
sumptions (ii) and (iii) polar oscillations decouple from axial ones 
in the linear regime. Therefore, we concentrate on purely axial os- 
cillations and evolve the (/3-component of the evolution variables 
only. We assume a conformally flat metric rl 



to the oscillations, which are related to the fluid 4-velocity accord- 



ds 



-a dt + (f> (dr + r d6 +r sin6 dip 



(1) 



where a is the lapse function and (f> the conformal factor, and con- 
sider a stress-energy tensor T''" of the form 

-'- ^ -* fluid + -^ mag + -* clas 



-2msS' 



(2) 



where p is the rest-mass density, h the specific enthalpy, P the 
isotropic fluid pressure, u'^ the 4-velocity of the fluid, b'^ the mag- 
netic field measured by a co-moving observer (with b^ :— b^b^), 
E'"' the shear tensor, and /is the shear modulus. The latter is ob- 
tained according to |Sotani et al.| ( (2007] l. 

The conservation of energy and momentum V^T*"^ = 0, and 
the induction equation lead to the following system of evolution 
equations 



= 0, 



(3) 



where g and 7 are the determinants of the 4-metric and 3-metric, 
respectively. The two-component state and flux vectors are given 
by 
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F" 



[S^, B" 
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w 



2^J.s^\ 



-v'^B' 



fryB 



2msS%, -^v^B 



(4) 
(5) 

(6) 



where B' are the magnetic field components as measured by an Eu- 



lerian observer I Anton et al.|2006f , and W — au^ is the Lorentz 



factor. The shear tensor S"'' = l/2g"^'^- contains the spatial 
derivatives (denoted by a comma) of the fluid displacement ^"^ due 

^ This provides a very good approximation as our neutron star models are 
almost perfectly spherically symmetric except for very small deviations due 
to the presence of an axisymmetric magnetic field. 



ing to ^ 



uf I'u!' , where v"^ is the (/j-component of the 



fluid 3-velocity. Hence, the evolution of the spatial derivatives ^"^^ 
and if ^ is given by 



(rj,t - (a^'"),* = with fcG{r,e}, 



(7) 



We also need to provide boundary conditions. At the surface 
the radial derivative of the displacement has to vanish (^"^^ = 0), 
as we assume a continuous traction and vanishing surface currents. 
At the crust-core interface we demand the continuity of the par- 
allel electric field, which implies a continuous displacement if . 
Together with the continuity of the traction the latter leads to a re- 
lation between the radial derivatives of the displacement in the core 
and in the crust: CS>ro,r = (1 + (J) C^rust.r with 5 = p.s/{brb''). 

To construct equilibrium models we choose between differ- 
ent barotropic EOSs for the core that are matched to an EOS for 
the crust. The available models for the core are the soft EOS A 
( Pan dharipandej 1 97 1 [ ), the intermediate EOS WFF3 l |Wiringa et al 
|1988] ), EOS APR jAkmal et al.|1998| l and the stiff EOS L ( |Pand' 
|haripande & Smith|1975| l. For the low density region of the crust 
we choose between EOS NV ([N egele & Vautherin 1973 ) and EOS 
DH (Douchm & Haensel|200lt . Details about the different combi- 
nations of EOSs can be found in |Sotani et al.lpOOT] ). We use a ref- 
erence model with a mass of IAA-Iq and a circumferential radius 
7?star = 12.26 km, described by the APR + DH EOS. In contrast 
to Sotani et al. (20071 we compute magnetised equilibrium models 
using the LORENE library (www . lorene . obspm. f r). 

Our simulation code is an extended version of the GRMHD 
code presented in [Cerda-Duran et al.H2009| . It includes the shear 
terms as they appear in (12), and the evolution of the displacements 
(jTl. The proper working of the MHD part of the code was demon- 
strated in [Cerda-Duran et al.H2008| l. To test the extended code we 
compared its results obtained for two limiting cases, zero magnetic 
field and zero shear modulus, with those of previous studies. The 
purely crustal shear oscillations presented in [Sotani et al.[ ( (2007] l are 
recovered with an agreement of 1%, and the Alfven continuum is 
obtained naturally as in [Cerda-Duran et al.[ ( |2009l l. Further details 
on the derivation of the model equations, the numerical methods, 
and the code tests will be discussed in Gabler, Cerda-Duran, Font 
& Stergioulas (in preparation). 

3 DAMPING OF CRUSTAL SHEAR MODES 

To study the behaviour of coupled crust-core oscillations we per- 
turb the equilibrium stellar model by imposing a velocity pertur- 
bation and then follow the time evolution of the system pi - (JTJl. 
Unless stated otherwise, we use a grid of 150x 100 zones in our 
simulations covering a domain [0, /Jstar] x [0, tt]. The angular grid 
is equidistant, while the radial grid is equidistant only in the crust, 
where 40 per cent of the zones are located, and coarsens towards 
the centre of the star. Symmetries are exploited whenever a pertur- 
bation is of purely odd or even parity with respect to the equatorial 
plane. We use the term damping in the following to refer to resonant 
absorption of crustal shear oscillations by the Alfven continuum of 
the core (unless we explicitly refer to numerical damping caused 
by finite-differencing). 

We investigated perturbations of different radial extent encom- 
passing only the crust, only the core, or the whole star. Since all 
three types of perturbations give qualitatively similar results, we 
focus on whole star perturbations in the following, as this is the 
most generic case. We first consider initial perturbations consist- 
ing of a single torsional, spherical vector harmonics Z-mode. To 
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Figure 1. Time evolution of overlap integrals with the eigenmodes of the crust. Left panel: Damping oil = 2,3, and 9 initial perturbations due to resonant 
absorption of the fundamental (n = 0) crustal shear mode for a magnetised model with 5 X 10^"' G (dots). In the corresponding unmagnetised models (solid 
lines) only numerical damping occurs which increases with the angular order I of the mode. Middle panel: Overlap integrals for the 1 = 2 mode, where 
resonant absorption of the crustal modes becomes stronger with increasing magnetic field strength. Right panel: same as middle panel, but for I = 9. 
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r [s] for B = 


5.500 


4.620 


0.260 


0.170 


T [s] for B = 5 X lO^-' G 


0.072 


0.080 


0.087 


0.081 


r [s] for B = IQi" G 


0.040 


0.040 


0.022 


0.021 



Table 1. Damping timescales r due to resonant absorption of crustal shear 
modes by the Alfven continuum for various initial perturbation modes I. 



investigate the damping of a single crustal mode we compute over- 
lap integrals of the evolved variables with mode eigenfunctions 
(the latter are found by solving the linear eigenvalue problem for 
the crustal modes, see.Schumaker & Thorne ( 1983 1; Messios et al.| 
( |2001^ ; [Sotani et ar] l |2007^ ). Because the eigenmodes of the crust 
form a complete orthonormal set we can expand any perturbation 
in terms of the corresponding eigenfunctions. The expansion fac- 
tors, which provide a measure of how strong each crustal mode 
contributes to the perturbation, are obtained via the overlap inte- 
grals with the eigenfunctions. For more details on this method see 
[Gabler et al.| ( |2009) . In the left panel of Fig.fTlwe show the max- 
imum (absolute) amplitudes of the overlap integrals for different 
initial perturbations and for simulations both without magnetic field 
(solid lines) and with a polar magnetic field of 5 x 10^'^G (dots) . 
In the field-free case the lines represent the numerical damping of 
crustal modes due to finite-differencing. When a magnetic field is 
present, the damping (now due to resonant absorption) increases 
with the magnetic field strength. 

For all modes, the timescale of resonant absorption is much 
shorter than that of numerical damping (see Table fTl. After about 
500 ms, the overlap integrals no longer sample the crust oscilla- 
tions, but instead the magneto-elastic oscillations which then dom- 
inate the evolution (see below). 

The middle and right panels of Fig.fT] show the overlap inte- 
grals for I = 2 and I = 9 modes of the crust as a function of time 
for different magnetic field strengths, respectively. For I = 2 and 
B = 5 X 10^^ G we find almost complete damping of the crustal 
mode after ~ 0.5 s. For a stronger magnetic field (B = 5 x 10^* G) 
the crustal mode becomes already damped after less than one oscil- 
lation, and only the dominant magneto-elastic oscillations remain. 
Similar statements hold for the evolution of the I = 9 perturbation. 
However, in that case it takes a few oscillations before the crustal 
mode is damped, even for B = 5 x 10^^ G. 

We also analysed a more general initial perturbation consist- 
ing of a mixture of Z = 2 up to Z = 10 modes, which excites a 
large number of crustal modes of different angular order I and ra- 
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Table 2. Damping timescales r due to resonant absorption of crustal shear 
modes by the Alfven continuum for initial perturbation modes I = 2,1 = 3 
and I = 9 and for different combinations of equations of state at B = 
10^* G. The number in the labelling of the EOS represents the mass of the 
neutron star model in Mq 
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Figure 2. Evolution of av'^ in the crust at 9 = 7r/4 for initial data contain- 
ing a large number of different perturbation modes. For vanishing magnetic 
field, only crustal modes are excited. For B = 5 X 10^"' G the crustal os- 
cillations are strongly damped, and for a ten times stronger magnetic field 
(B = 5x10^* G) magneto-elastic oscillations dominate the evolution from 
the start. 



dial order n. Fig.[2] shows the resulting evolution of the velocity in 
the crust at ^ = n/4. For the unmagnetised model, low frequency, 
fundamental (n = 0) oscillations can easily be distinguished from 
higher frequency overtones with n ^ 1. When increasing the mag- 
netic field strength to 5 x 10^"^ G the lower frequency modes are 
completely damped after ~ 250 ms, whereas the higher frequency 
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Figure 3. Spatial distribution of Fourier amplitudes at several QPO frequencies. Three different types of QPOs are present: (i) upper QPOs in panels a, d, f, g, 
i, and k; (ii) lower QPOs in panels b and e; (iii) edge QPOs in panels c, h, j, and /. The two black dashed lines mark the location of the crust, while the blue 
lines represent magnetic field lines. The colour scale ranges from zero amplitude (white) to maximum amplitude (black). 



overtones survive for a longer time. In the long run, a low frequency 
(~ 0.3 Hz) magneto-elastic oscillation dominates. At the largest 
magnetic field strength shown here, 5 x 10^* G, there is no sign 
of either low or high frequency crustal modes, the evolution being 
completely dominated by magneto-elastic oscillations. 

Damping timescales of n = 0, / = {2, 3, 9} modes for dif- 
ferent EOSs are listed in Table [2] For all EOSs and masses stud- 
ied here, we find significant damping of the crustal shear modes at 
B = 10^* G. Up to the highest-order / — 9 mode, the damping 
timescales are shorter than 0.1s. From the results in Tables 1 and 
2, we deduce that even at B = 5 x lO^"' G, which is up to two or- 
ders of magnitude weaker than assumed magnetar field strengths, 
the damping timescale is shorter than 0.2s, which is more than two 
orders of magnitude shorter than the duration of observed QPOs. 

4 LONG-TERM QPOS 

Besides the damping of crustal modes, we observe long-lasting os- 
cillations in the fluid core of the magnetar. These long-term QPOs 
are identified by local maxima in Fourier space. Let us consider 
an intermediate magnetic field strength of 4 x 10^* G, where both 
the magnetic field and the crust influence the dynamics. As in the 
case without crust ( |Sotani et al.|2008b[ |Cerda-Duran et al.|2009) 
we find two different families of long-term QPOs, as demonstrated 
in Fig.lS] which shows the spatial distribution of Fourier amplitudes 
at several QPO frequencies. The lower QPOs (LJi ) are located in- 
side the region of closed field lines, while the upper QPOs (U„ ) 
concentrate along open magnetic field lines closer to the poles. We 
computed QPOs of either odd (— ) or even (+) parity w.r.t. the 
equatorial plane, which allows for a better identification of QPOs 
of similar frequency but opposite parity. 

The lower QPOs (Fig.lSl panels b and e) appear to be similar to 
those found for models without crust, except that they are limited to 



the region of closed magnetic field lines inside the core. The upper 
QPOs are influenced by the presence of the crust in several ways. 
First, they are limited to the fluid region, and become vanishingly 
small at the base of the crust (Fig.lS] panels a, d,f, g, i, and k), where 
they are reflected. This behaviour is similar to that caused by the 
boundary conditions in jSotani et al.| ( (2008b^ , and differs from that 
of the pure fluid case considered in |Cerda-Duran et al.| ( |2009] ). In 
the latter work the continuous traction boundary condition imposed 
at the surface of the star resulted in a strong displacement there. 
Their different behaviour at the base of the crust (w.r.t. the pure 
fluid case in |Cerda-Diiran et al.| l |2009l ) causes a rearrangement of 
the QPOs. Here, the lowest frequency QPO (panel a) is symmetric 
w.r.t. the equatorial plane (even parity) - while it did not exist at all 
in |Cerda-Duran et al.| ( |2009j - and the lowest-frequency QPO has 
odd parity. 



While QPOs are located close to the symmetry axis of the field 
(polar axis) in models without crust ( Sotani et al. 2008b[ |Cerda-| 
|Duran et al.||2009| ), they are attached to field lines crossing the 
equator at around 4 km in our models including the crust. Appar- 
ently the strong coupling introduced by the crust complicates the 
oscillatory behaviour of the field lines, such that the interaction be- 
tween neighbouring polar field lines prevents the QPOs from being 
established. 



Furthermore, we find a new family of QPOs (Fig.lS] panels c, 
h, j, and Z) connected to the last open field line of the fluid core, 
each member representing the lower-frequency edge of an Alfven 
continuum along the open field lines. QPOs at similar locations 
were also identified in simulations by |Colaiuda et al.| ( |2009| ) without 
a crust. 
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5 DISCUSSION 

In this Letter we have presented the first numerical simulations of 
axisymmetric, torsional, magneto-elastic oscillations in a realistic 
magnetar model, including an extended crust. 

We focus on the timescales for resonant absorption of low 
frequency, n — crustal shear modes by the Alfven continuum 
of the core and find that even at B = 5 x 10^"^ G, which is up 
to two orders of magnitude weaker than assumed magnetar field 
strengths, the damping timescale is shorter than 0.2s, which is more 
than two orders of magnitude shorter than the duration of observed 
QPOs. Furthermore, the crust EOSs NV and DH which we have 
used have a very high shear modulus, compared to a range of other 
proposed EOSs ( [Steiner & Watts|2009^ . Comparing the damping 
timescales obtained here for the NV and DH EOSs, we find that, 
at a given magnetic field strength, a lower shear modulus results 
in a shorter damping timescale due to resonant absorption. Thus, 
for lower values of the shear modulus than considered here, the de- 
duced damping timescales would be even shorter than our current 
findings. In addition, no significant excitation of crustal modes by 
the Alfven continuum was observed at the end of our simulations, 
which reached up to several seconds. The above results do not leave 
much room for shear modes of the crust to be able to sustain oscil- 
lations that lead to significant modulations in the X-ray tail of giant 
SGR bursts lasting for several tens of seconds, when the magnetic 
field is global. The case of magnetic field configurations confined 
to the crust, which may be realised if the core is a type I supercon- 
ductor, requires further investigation. 

In contrast to the shear modes, the Alfven QPO model (see 
Levi n (2007); Sotani et al. ( 2008b) ; [Cerda-Duran et a"L] ( |2009l l; |Con 
laiuda et al.|(|2009|l) has several attractive features, matching to ob- 



served frequencies for a dipole field strength up to several times 
10^^ G and explaining the observed integer ratios for the frequen- 
cies of QPOs. Here we find that for magnetic field strengths less 
than about 10^^ G, the Alfven QPOs are mostly confined to the 
fluid core. For EOSs with a smaller shear modulus, this could re- 
duce to several times 10^* G. Above this minimum field strength, 
Alfven oscillations completely dominate in the crust and the result- 
ing QPOs have a maximum amplitude at the surface of the star. 
It is thus likely, that the Alfven QPO model can operate efficiently 
only above such a minimum field strength. In that case, the possible 
strength of magnetar magnetic fields is limited to a narrow range, 
between the minimum strength discussed here and the upper bound 
determined in |Cerda-Duran et al.| ( [2009l >. Furthermore, the empir- 
ical formulas for the QPO frequencies presented in |Cerda-Duran| 
|et al.| ( [200 9^ could be used to directly constrain the strength of the 
magnetic field and the compactness of the star. 

A more extended investigation of the minimum magnetic field 
strength for the Alfven QPO model will appear in |Gabler et al.| 
(J20T0l. We are also planning to extend our model by taking into 
account additional effects, such as different magnetic field topolo- 
gies (see also Sotani et al.| pOOSal), the coupling of the interior 
dynamics to an external magnetosphere, and the effect of field-line 
entanglement in the core. 
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